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ABSTRACT 

Changes  in  round-trip  acoustic  travel  time  (r)  measured  between  a  bottom  moored  inverted  echo  sounder 
and  the  sea  surface  can  be  interpreted  as  changes  in  dynamic  height  (/>)  with  suitable  calibration  information. 

The  t.  D.  and  isotherm  and  isopycnal  depths  (Z)  have  been  calculated  using  hydrographic  (CTD)  data  from 
three  regions:  the  Norwegian  Current,  the  Sargasso  Sea  and  the  eastern  Gulf  Stream.  Regressions  of  l)  and  / 
on  J  were  performed.  The  slope  for  the  Norwegian  Current  is  -  1 .70  +  0.01  [dyn  cm  (m  s)~t],  for  the  Sargasso 
Sea  -3.89  +  0.16  [dyn  cm  (m  s)’T  and  for  the  Gulf  Stream  -3.13  ±0.07  [dyn  cm  (m  sL'J.  The  quasiirandom 
scatter  about  regression  curves  is  found  to  be  primarily  the  result  of  variability  in  the  seasonal  thermocline. 
where  present. 


1.  Introduction 

The  inverted  echo  sounder  (IES:  Bitterman  and 
Watts,  1979)  is  a  bottom  moored,  upward-looking 
acoustic  device  for  measuring  the  round-trip  acoustic 
travel  time  (t)  between  the  sea  surface  and  bottom. 
Fluctuations  in  the  height  of  the  sea  surface  have  a 
measurable  effect  on  r  ( 1  ms  m  " Watts  and  Rossby. 
19771  but  the  dominant  factor  by  far  is  variability  of 
the  second-speed  profile  in  the  intervening  water  col¬ 
umn.  For  example,  a  change  in  the  thermal  structure 
corresponding  to  a  change  in  dynamic  height  (D)  of 
100  dynamic  centimeters  (dcm)  results  in  a  change  in 
r  of  about  30  ms  in  the  Gulf  Stream  region  (equivalent 
to  about  a  600  m  change  in  thermocline  depth).  The 
IES  is  thus  well  suited  for  the  study  of  mesoscale  baro- 
clinic  variability.  If  the  IES  is  also  provided  with  a  pre¬ 
cision  pressure  sensor  (Watts  and  Wimbush.  1981),  the 
combined  measurements  yield  both  barotropic  and 
baroclinic  components  of  variability  in  the  surface 
pressure  field. 

Interpreting  changes  in  IES  measured  r  as  changes 
in  D  (or  pycnocline  depth)  requires  supporting  infor¬ 
mation  about  the  sound-speed  and  density  profiles  in 
the  region  where  the  IES  is  to  be  deployed.  This  infor¬ 
mation  can  be  obtained  from  hydrographic  observa¬ 
tions  at  IES  sites  wherever  practicable  and/or  from  sta¬ 
tistics  of  historical  hydrographic  data  from  the  region 
of  interest.  The  study  described  in  this  report  is  an 
example  of  the  latter. 

Calculations  of  r,  D ,  and  isotherm  and  isopycnal 
depths  were  made  using  hydrographic  (CTD)  profile 
data  acquired  in  three  regions  where  mesoscale  features 
were  present:  the  southeastern  Norwegian  Sea  (NS), 
the  eastern  Gulf  Stream  (GS)  south  of  New  England, 
and  the  Sargasso  Sea  (SS)  east  of  Jacksonville,  Florida. 


The  results  are  based  exclusively  on  these  data  (no  ac¬ 
tual  IES  data  are  included). 

2.  Background 

The  attempt  to  use  bottom-mounted  inverted  echo 
sounders  (IESs)  to  monitor  changes  in  the  depth  of  the 
main  thermocline  was  first  made  by  Rossby  (1969). 
Using  two  hydrographic  datasets,  he  compared  cal¬ 
culated  round-trip  acoustic  travel  lime  (r)  with  the 
depths  of  the  10°  and  I2°C  isotherms  and  found  re¬ 
markable  correlations.  A  linear  least-squares  slope  of 
about  20  m  (ms)  1  was  found.  He  subsequently  con¬ 
ducted  an  eight-day  experiment  near  Bermuda,  where 
he  found  that  measured  r  was  quite  stable,  varying 
only  in  response  to  surface  and  internal  tides.  During 
MODE,  several  bottom-moored  IESs  were  deployed 
to  monitor  thermocline  depth  fluctuations.  There  were 
simultaneous  hydrographic  (CTD)  observations  in  the 
area.  Watts  and  Rossby  (1977)  summarized  the  MODE 
results  and  formally  examined  the  relationship  between 
r  and  D.  They  found  that  fluctuations  in  D  could  be 
inferred  from  t  changes  with  an  uncertainty  of  only 
±  1  dcm.  To  arrive  at  this  estimate,  they  intercompared 
r  from  the  IES  with  r  and  D  calculated  from  CTD 
data.  These  -alculations  were  done  across  the  main 
thermocline  between  500  and  1 500  m.  which  excluded 
the  surface  layers.  Comparisons  of  measured  r  with 
the  depth  of  the  10°C  isotherm  (thermocline  depth) 
showed  a  rms  scatter  of  about  5  m,  which  scales  to 
about  0.7  dcm  in  D.  Thermocline  depth  varies  as  /), 
provided  the  first  baroclinic  mode  dominates  the  vari¬ 
ability.  They  also  calculated  t  and  D  between  the  sur¬ 
face  and  4000  m.  using  historical  hydrographic 
data  from  the  Gulf  Stream  region,  and  obtained  an 
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rms  uncertainty  of  5  dyn  cm  about  a  linear,  least- 
squares  fit. 

3.  CTD  data  analysis 

a.  Approach 

Changes  in  D  and  r  are  calculated  from  profiles  of 
specific  volume  («),  and  reciprocal  sound  speed  (7), 
which  are  functions  of  pressure  (/’),  temperature  ( T ) 
and  salinity  (5).  Integrating  in  pressure  coordinates 

D'  =  -['a'dP  (1) 

Jr, 

2  fF- 

t'~  —  I  ( ay)'dP ,  (2) 

gJr, 

wnere  D  =  D  +  D\  r  =  f  A  r':  the  overbar  denotes  a 
local  time  average. 

The  integrands  of  ( I )  and  (2)  are  of  similar  form  and 
the  arguments  presented  by  Watts  and  Rossby  (1977) 
demonstrate  how  D'  and  r'  are  nearly  linearly  related. 
The  relationship  between  D' and  t' depends  on  hydro- 
graphic  conditions  in  the  vicinity  of  an  1ES  mooring, 
and  ideally,  we  should  like  to  have  a  series  of  hydro- 
graphic  profiles  at  each  IES  site  to  provide  stable  sta¬ 
tistics.  In  a  region  such  as  the  Gulf  Stream,  which  varies 
on  time  scales  ranging  from  several  days  to  several 
months,  this  is  not  practical.  In  such  a  region,  however, 
variability  is  primarily  due  to  lateral  shifting  of  hori¬ 
zontal  gradients  and  the  necessary  statistics  can  be  ob¬ 
tained  by  quasi-sy  noptic  sampling  over  the  region  sur- 
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Flo.  2.  As  in  Fig.  I  but  between  100  and  1  (XX)  db 


rounding  the  IESs.  There  are  seasonal  warming  and 
cooling  effects  which  may  be  missed  with  this  approach, 
but  the  svnoptic-scale  data  can  be  complemented  by 
representative  XBT  and/or  AXBT  surveys  throughout 
the  year. 

In  the  analysis  to  follow.  CTD  data  from  three  re¬ 
gions  are  examined.  D',  t'.  and.  in  some  cases,  isopleth 
depths  are  calculated  and  linear  regressions  are  per¬ 
formed.  Sources  of  uncertainties  and  scatter  in  the 
regressions  are  discussed. 


h.  Description  of  the  data 


Fki.  I  P  vs  1 
calculated  with  (’ 
20  and  1000  db 


"  (uM 

.  with  superimposed  quadratic  regression  curve. 
TD  profiles  from  the  Norwegian  Current,  between 


Area  1  is  the  southeastern  Norwegian  Sea  (NS),  and 
the  data  consist  of  48  CTD  casts  to  about  1 1 00  m, 
taken  in  three  sections  across  the  Norwegian  Current. 
This  current  is  the  extension  of  the  North  Atlantic 
Current  and  transports  Atlantic  water  northward 
through  the  Faeroe-Shetland  Channel  to  mix  with  wa¬ 
ters  of  Arctic  and  Subarctic  origin.  The  data  were  ac¬ 
quired  in  March  1981  when  no  seasonal  thermocline 
was  observed.  The  velocity  of  the  Norwegian  Current 
was  about  25  cm  s  1  from  the  surface  to  a  depth  of 
about  300  m.  (These  data  were  obtained  by  the  Naval 
Oceanographic  Office  (NOO)  and  have  been  part  of 
NORDA's  archives.) 

Area  2  is  the  western  Sargasso  Sea  (SS)  near  30°N. 
7 1 0  W.  Data  consist  of  a  grid  of  about  30  CTDs  to  2000 
m.  separated  by  about  20  nautical  miles  (n  mi),  ac¬ 
quired  in  September  1979.  The  area  contained  an  an¬ 
ticyclone;  eddy  with  velocities  of  the  order  of  20  cm 
s  1  ( Hallock  et  al..  1981)  and  there  was  a  well-developed 
seasonal  thermocline  near  60  m.  The  source  of  these 
data  is  the  same  as  for  Area  1 . 
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f'H  t.  }  Vs  in  Fig.  1  bui  tjlculjifj  wnb  12  ('l  l)  profiles 
from  the  Sargasso  Sea  between  12  and  1^50  db 


I  K i  5  \s  in  Fig  1  bin  calculated  with  <><>  ('ll)  profiles 
from  the  (»u!f  Stream  between  *  and  VMKi  Jb 


Area  3  consists  of  6(.  CTD  casts  from  the  eastern 
Gulf  Stream  region  (OS )  in  the  vicinity  of  the  New 
England  Seamounts.  The  data  are  distributed  in  sec¬ 
tions  across  and  along  the  Gulf  Stream  and  most  extend 
to  the  bottom.  A  seasonal  thermocline  was  present  at 
most  locations  in  and  south  of  the  Stream.  There  were 
strong  frontal  features  near  the  northern  boundary  of 
the  current.  These  data  were  acquired  by  the  Woixis 
Hole  Oceanographic  Institution  in  summer  1983.  and 
are  described  by  Hogg  et  al.  ( 1985). 


c  Cah'uUinom 

I  or  each  of  the  three  datasets.  /)  and  -  wcie  com¬ 
puted  by  integrating  Tqs  1 1 )  and  (2)  oxer  selected  pres¬ 
sure  intervals.  The  pressure  limits  (/'  and  /’  I  were 
chosen  to  include  the  maximum  tuimbet  of  proltles. 
and  >ct  get  close  to  the  surface  as  ss oil  as  below  most 
of  the  baroclinic  xaruhiltty .  A  second  set  o!  calculations 
was  done  ot  l>  and  .-  where  the  upper  pressure  limit 
was  chosen  deliberately  to  exclude  the  upper  layer  and 
seasonal  thermocline.  where  present.  />  is  plotted  as  a 


Fto.  4.  As  in  Fig.  3  but  between  200  and  1950  db. 
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Fni  6  As  in  Fig  5  but  between  200  and  3000  db 
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function  of  r'  in  Tigs.  1-6  whore  quadratic  regression 
curves  have  been  superimposed 

1  inear  (,)  and  quadratic  (,()  regressions  of/)  and  r' 
and  error  estimate  (/•.)  were  computed  as  follows 

/)  -  K  t ' ■+  ( '  (linear) 

t.  rms(/>  /)  I 

/)„  I,.t  '  -  H.t  *  (  .,  (quadratic) 

/  rmsi/)  /)„. ) 

I  he  uncertainty  in  the  linear  slope.  H  .  is  given  hv 
H  /.',/[  V1  rmsi  -  )]  where  N  is  the  numher  <>t  obser¬ 
vations. 

Results  of  the  regression  calculations  for  />  vs  r  are 
summarized  in  I  ables  t  and  2  (since  />'  and  •  are 
anomalies.  (  -  0  and  is  not  shown)  I  mcar  regressions 
were  also  done  for  selected  isotherm  and  isopvcnal 
depths  vs  r  and  these  results  are  shown  m  fable  .3 

(J  DiscusMon 

It  is  apparent  in  all  cases  that  a  sigmheant  correlation 
exists  between  D  and  r  .  Indeed,  the  values  of  rare  all 
greater  than  0.97.  The  departures  from  a  perfect  linear 
lit  fall  into  two  categories:  a  general  curvature  which 
can  be  corrected  for  with  a  quadratic  fit.  and  a  quasi¬ 
random  residual  scatter. 

The  NS  area  data  yields  the  best  fit  of  P  and  r‘.  and 
case  NS-2  shows  very  little  improvement  over  NS-1. 
A  significant  seasonal  thermocline  is  absent  and  near¬ 
surface  variability  is  that  of  salinity  in  a  thin  (10  m) 
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I  he  slope.  H  . 

is  about  half  the 

values  found  in 

the  other  two  regions  due  to  a  significantly  different 
/-.S'  relationship 

I  he  SS  data  exhibits  the  lowest  value  of  r.  I  he  sari 
abilitv  of/)  is  comparable  to  that  in  NS.  but  the  signal- 
to-noisc  ratio  is  lower  There  is  a  well-developed  sea 
sonal  thermocline.  as  well  as  a  subsurface  maximum 
of  density  variability  near  40  m  The  value  of  /  is 
slightly  greater  than  that  given  by  Watts  and  Rossbv 
( Id-1^)  for  the  same  region  1  xclusion  of  the  upper  2(Mi 
m  ((  use  SS-2  i  reduces  /■.  and  improves  the  lit  dra¬ 
matically 

Signal  levels  in  the  ( iS  area  data  are  about  eightfold 
greater  than  m  the  othei  areas  I  he  /  is  larger  than 
lor  SS  but  •  is  considerably  higher,  indicating  a  highvi 
signal-to-noise  ratio  \  seasonal  thermocline  is  present 
in  and  south  ol  the  (  mil  Stream  and  a  subsurface  mav- 
imuni  in  density  sai  lability  extends  from  about  V'  to 
65  m  Case  GS-2  reduces  /  bv  about  me  half  and 
improves  the  tit  significantly 

The  results  of  the  quadratic  regressions  show  a  gen¬ 
eral  improvement  for  all  cases.  This  improvement  is 
due  to  the  accommodation  of  the  small  but  systematic 
curvature  in  /)'  vs  r  '.  The  curvature  effect  is  estimated 
by  . t^r '■  (which  is  numerically  equal  to  ( since  r 
and  /)'  have  zero  mean)  and  is  of  the  same  order  as 
the  rms  error  of  the  quadratic  fit.  A  statistical  test 
(/•-Test)  shows  that  significant  improvement  results 
with  the  quadratic  term  for  all  six  cases  The  impiove- 
ment  is  marginal,  however,  for  SS-2.  In  the  NS  cases, 
the  curvature  effect  is  nearly  negligible.  In  both  the  GS 
and  SS  cases  however,  the  curvature  effect  as  well  as 
/■..'„  are  considerably  reduced  by  calculating  D'  and 
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below  the  thermocline.  but  the  slopes  Hi.  Hu  do  not 
change  significantly  The  source  of  the  curvature,  which 
does  not  always  have  the  same  sign,  is  a  result  of  the 
nonlinearities  of  both  the  equation  of  state  and  the 
I  -S  relationship. 

In  cases  where  the  first  baroclmic  mode  dominates 
variability,  dynamic  height  changes  are  reflected  in 
isotherm  and  isopycnal  depth  changes  in  the  main 
thermocline.  When  IFS  data  are  used  for  monitoring 
positions  of  strong  mesoscale  features  or  w  here  a  two- 
layer  description  is  sought,  a  direct  interpretation  of  r' 
as  changes  in  thermocline  depth  (/,)  or  pycnocline 
depth  (/,)  may  be  useful  C  orrelations  between  /;. 

and  r  (  Table  3:  the  tsopycnals  selected  corresp  aid 
approximately  to  the  isotherms  indicated  by  the  tem¬ 
perature  in  parentheses)  are  generally  lower  than  those 
lor  /)  and  r  .  This  is  partly  because  //.  /.,  are  likely 
to  he  noisier  quantities  than  /T(an  integral)  and  partly 
because  events  at  some  other  depth  may  affect  t '  but 
not  a  particular  isotherm  (or  isopycnal). 

Tor  the  SS  case.  «/(/ 1,,>  and  /  .(/,,.)  are  24  m  (ms)  ’ 
and  10  m.  respectively,  which  differ  considerably  from 
the  14  m  (ms)  1  and  4  m  reported  by  Watts  (1975). 
for  the  same  region.  This  difference  might  be  seasonal 
or  regional  (Watts'  data  was  acquired  in  summer  near 
28°N.  69.5°W  while  the  SS  data  was  acquired  in  early 
fall,  near  30°N.  7 1  “W’l.  For  both  CiS  and  SS.  /,  above 
the  main  thermocline  is  poorly  correlated  with  r':  this 
observation  may  have  implications  for  relations  be¬ 
tween  near-surface  and  mesoscale  (deep)  variability 
Results  for  /,  in  the  Norwegian  C  urrent  are  nearly 
the  same  for  all  isotherms  The  vertical  temperature 
gradient  is  nearly  constant  between  the  surface  and  the 
lower  depth  limit  of  the  data,  and  salinity  variability 
is  low.  In  all  cases.  /...  t'  has  a  slightly  lower  correlation 
than  does  /,.  r  at  the  corresponding  isotherm  This 
is  expected,  since  sound-spiced  fluctuations  tend  to  fol¬ 
low  temperature  fluctuations  more  closely  than  do 
density  fluctuations. 

The  /.(NS).  /i;  (SS)  and  /,„ (CiS)  have  the  highest 
correlations  with  r1  and  arc  the  most  useful  for  inter¬ 


preting  1ES  observations.  Also  /i;  and  /,■.  have  good 
correlations  and  are  included  since  they  have  been  used 
to  define  the  thermocline  by  other  investigators  le  g 
Hansen.  1964:  Watts  and  Johns.  1982;  Watts  and 
Tracey.  1985). 

To  determine  possible  effects  of  the  lower  pressure 
limit.  /’, .  calculations  with  a  subset  of  the  CiS  data 
which  extended  to  5000  db  were  performed  between 
12  and  3000  db  and  between  12  and  5(XM)  db.  I  here 
is  very  little  difference  between  the  results  in  these-  two 
cases 

The  salient  results  to  emerge  from  these  calculations 
are  the  values  of  H,  and  their  associated  errors  1  hex 
are  the  regional  calibration  factors  for  conversion  ot 
II  S  measured  r'  into  /) .  / ,  or  1  he  primarv  lim¬ 
itation  to  the  accuracy  of  the  vs  r  regression  tor 
any  given  area  appears  to  be  variability  in  the  upper 
100-200  m  of  the  water  column. 

c  further  analysis  ot  saltier 

In  cases  where  a  seasonal  thermocline  is  present,  the 
upper  200  m  is  i  major  source  of  scatter  in  the  />  vs 
r  relation  It  is  useful  to  explore  further  the  nature  ot 
this  scatter  I  he  NS  cases  contain  relatively  little  scatter 
and  hence  can  serve  as  a  control  tor  comparison  with 
the  other  two  areas 

•\  linear  relationship  between  />  and  r  requires  some 
similarity  in  the  integrands  of  ( 1 )  and  (2)  \  necessarv 
(but  not  sufficient)  condition  is  that  profiles  ot  rmslu  > 
and  rms  (toy  )'|  have  the  same  shape  I  hese  profiles  for 
the  three  regions  are  plotted  in  Tigs  7-4  1  here  is  a 
subsurface  maximum  mu'  for  both  ( iS  and  SS.  w  Inch 
is  not  present  in  (<«>V  1  his  feature,  which  corresponds 
to  the  seasonal  thermocline.  gives  rise  to  much  of  the 
scatter  observed  in  the  />  -r  relationship  1  his  is  due 
to  a  different  IS  relationship  in  the  upper  lasers 
where  salinity  variations  tend  to  change  density  in  the 
same  sense  as  temperature  variations,  i.e  .  7  and  .5  are 
noncompensating.  No  subsurface  maximum  in  o  is 
present  in  NS;  a  thin  layer  of  low  salinity  near  the 
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surface,  a  manifestation  of  the  fresh  Norwegian  coastal 
water,  has  xerx  little  effect  on  the  integral  of  <*'.  The 
shallow  xanabilitx  in  o  does  not  appear  to  be  correlated 
with  the  primary  signal  which  corresponds  to  the  xari- 
abilitv  in  the  main  thermochne  t  his  shallow  variability 
is  likelx  to  be  spatiallv  and  temporalis  aliased:  it  mas 
be  associated  ssnh  internal  ssases  or  frontal  lluctua- 
Uons 

1 o  better  understand  the  effects  ol  temperature  and 
salinity  sanations  on  the  /)  -  relation,  sse  focus  at- 
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tention  on  CiS  data  and  examine  Iqs  (II  and  (2)  more 
closely  I  he  integrands  can  be  expanded  about  the 
mean  to  yield 
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t'  1  7,1  *  7/ :  7  *M  *  7s:  "  T !  *  rs. 

These  quantities  were  computed  for  the  GS  data  oxer 
three  pressure  intervals  (3-3000.  3-200.  200-3(XX)  db) 
and  their  statistics  appear  in  Table  4.  The  dominant 
terms  in  the  P'-t  relation  are  P) .  P\.  and  t;i.  Cor¬ 
relations  are  slightly  higher  among  all  terms  for  the 
200-3000  db  integral  than  for  3-3000  db.  They  are 
significantly  lower  for  3-200  db.  indicating  a  different 
regime  in  this  interval,  as  previously  suggested. 

In  the  deeper  interval  /)s  0.4 1 P, :  the  high  neg¬ 

ative  correlation  implies  the  effect  of  sail  nits  vanabilit> 
is  almost  exactly  opposed  to  that  of  temperature  vari¬ 
ability.  Also  in  this  interval.  t\  *  0.07r;  and  r 
-  O.OSr / , .  showing  the  dominance  of  rn.  In  the 
shallower  interval,  correlations  retain  their  signs  but 
are  significantly  lower.  /)s  •=.  0.52 P,.  rs  =  O.lbr, 

and  t,i  *  0. 1  3t  , , .  show  ing  the  increased  importance 

of  salinity  vanahility  to  r'.  In  all  cases.  rs,  rs:  with 
high  correlation  The  increased  effect  of  salinity  and 
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poorer  T-S  correlation  in  the  upper  200  m  indicate 
that  seasonal  corrections  based  only  on  temperature 
profiles  (e.g..  XBTs)  will  have  limitations  in  this  in¬ 
terval. 

4.  Summary 

The  GS-2  and  SS-2  cases  simulate  an  extreme  sea¬ 
sonal  variation,  since  choosing  PI  at  200  db  is  roughly 
equivalent  to  extending  the  T  and  S  at  200  db  to  the 
surface.  Hence,  seasonal  changes  in  Bt  of  the  order  of 
those  shown  in  Table  1  might  be  expected.  Clearly,  a 
study  of  year-round  observations  would  provide  more 
accurate  estimates  of  seasonal  variations. 

If  one  is  interested  primarily  in  the  dynamic  signal 
in  the  main  thermocline.  as  suggested  in  the  calcula¬ 
tions  done  by  Watts  (1975)  (see  section  2.  this  paper). 
H i  for  the  deeper  PI  might  be  used,  with  its  better  sta¬ 
tistics.  thus  avoiding  the  seasonal  variation.  Further¬ 
more.  additional  regression  calculations  with  SS  and 
GS  data  (not  presented  here)  with  PI  at  200  for  D'. 
but  with  PI  remaining  near  the  surface  for  r'.  corre¬ 
lations  are  nearly  as  good  (at  least  0.99 1 )  as  for  SS-2 
and  GS-2.  On  the  other  hand,  if  the  main  interest  is 
specifically  in  the  dynamic  height  changes  at  the  sur¬ 
face.  some  estimate  of  seasonal  variations  should  be 
included  when  interpreting  r\ 

The  most  direct  dynamical  interpretation  of  IES 
travel-time  data  is  as  changes  in  dynamic  height  across 
the  main  thermocline,  when  mesoscale  variability  is 
present.  The  relationship  is  virtually  linear  and  uncer¬ 
tainties  are  acceptably  small  in  the  three  regions  ana¬ 
lyzed  here.  Comprehensive  hydrographic  surveys  can 
provide  the  necessary  statistics  to  effect  the  conversion 
of  travel-time  data  to  dyr  mic  height  change.  These 


surveys  may  be  concurrent  with  an  1F.S  deployment 
or  composed  of  historical  data  They  must  be  extensive 
enough  to  include  all  sets  of  hydrographic  conditions 
that  are  likely  to  be  experienced  at  the  IES  sites.  The 
slopes  and  uncertainties  for  the  NS.  SS  and  GS  regions 
are  1 .70  ±  0.01  (dem  ms  ').  3.X9  ±  0. 16  (dem  ms  ') 

and  3.13  ±  0.07  (dem  ms  ').  respectively.  There  is  a 
small  but  significant  quadratic  dependence  of /)' on  r  '. 
In  applications  where  high  precision  is  required,  results 
such  as  those  in  Table  2  will  be  useful. 

The  inference  of  isotherm  depth  change  from  travel- 
time  data  depends  more  strongly  on  the  assumption 
that  the  first  baroclinic  mode  is  dominant,  which  in 
SS  and  GS  is  likely  to  be  a  good  assumption.  For  the 
isotherms  which  gave  the  best  correlations  with  r'.  the 
slopes  and  uncertainties  for  the  NS  (/0.  SS  (/,;).  GS 
(/hi)  are  -35.2  ±  0.6  (m  ms  ’).  -22.5  ±  0.8  m  ms  1 
and  -21.1  ±  0.2  m  ms  '.  respectively. 

The  primary  source  of  uncertainty  in  the  conversion 
factor  is  relatively  uncorrelated  T-S  variability  in  the 
near  surface  layers  when  a  seasonal  thermocline  is 
present.  This  variability  is  likely  to  be  associated  with 
local  warming,  frontal  movement  or  internal  waves, 
and  is  probably  aliased  in  the  data  sets  of  this  study. 
In  future  calibration  studies,  attention  should  be  fo¬ 
cused  on  the  shallow  effects  with  data,  which  are  able 
to  resolve  the  smaller  scales. 
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